Anisotropic inverse Compton scattering from 
the trans— relativistic to the ultra— relativistic 
regime and apphcation to the radio galaxies. 



o\ ■ 

ON ■ G. Bmnetti ^ 

Dip. di Astronomia, Univ. di Bologna, via Ranzani 1, 1-40127 Bologna, Italy. 
. Istituto di Radioastronomia - CNR, via Gohetti 101, 1-40129 Bologna, Italy. 

o ■ 
c<i : 

^ , Abstract 

^ : 

. The problem of the anisotropic Inverse Compton scattering between a monochro- 

I matic photon beam and relativistic electrons is revisited and formally solved without 
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0> . and with a population of electrons isotropically distributed, under the assumption 

I that the energy distribution of the relativistic particles follows a simple power law as 

' it is the case in many astrophysical applications. Both the Thomson approximation 
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1 Introduction 



Inverse Compton (IC) scattering is one of the most common radiative processes 
in astrophysics. The papers of Jones (1968) and Blumenthal & Gould (1970) 
provide the basic theory, while both exact expressions and useful approxi- 
mations for interaction rates, mean scattered energies, and dispersion about 
the mean energies are given by Coppi & Blandford (1990). The computations 
reported in these papers have been performed by supposing an isotropic dis- 
tribution of the incident photon momenta. However, there are astrophysical 
situations in which the photon momenta are anisotropically distributed so that 
the standard equations are not sufficient. 

From a theoretical standpoint, Baylis et al.(1967) showed that the Compton 
spectrum (emitted over the total solid angle) from a beam of ultra-relativistic 
electron population (iVe(7) oc 7^*^) in the Thomson approximation is inde- 
pendent on the incoming photon distribution being for all considered cases 
oc e"*-^^^^/^. The anisotropic inverse Compton (AIC) equations describing the 
scattering between a photon beam and isotropic relativistic electrons in the 
Thomson approximation were obtained by Bonometto et al.(1970) for a power 
law energy distribution of ultra-relativistic electrons. Aharonian & Atoyan 
(1981) obtained the general redistribution function describing the scattering 
probability into a given direction of a monochromatic incoming photon beam 
scattered by isotropic monochromatic relativistic electrons. Furthermore, by 
integrating the redistribution function over a population of ultra-relativistic 
electrons with a power law energy distribution (iVe(7) oc 7"*^) they obtained 
the analytical expression of the spectrum (in terms of hypergeometric series) 
and showed that it presents the same shape of that derived from an isotropic 
illumination (i.e. oc e-(''-^)/2 in the Thomson approximation and oc e ^ln{eeo) 
in the Klein-Nishina regime; eg being the energy of the seed photons). More 
recently, Nagirner & Poutanen (1993) have obtained the redistribution ma- 
trix for Compton scattering of polarized incoming photon beams and, by in- 
tegrating the redistribution function over a power law energy spectrum of 
ultra-relativistic electrons in the Thomson approximation, they showed the 
equivalence between their results and those of Bonometto et al.(1970). 

All the formulae of astrophysical interest reported in these papers are obtained 
by making use of the ultra-relativistic limit. As a consequence, they are not 
useful to the detailed study of astrophysical situations in which also trans- 
relativistic and mildly-relativistic electrons are involved. 

More recently Brunetti et al.(1997) have treated the AIC Thomson scattering 
of a photon beam by an isotropic electron distribution not restricted to the 
ultra-relativistic limit. 
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The AIC has been considered in several specific astrophysical situations. The 
Compton losses suffered by electrons moving in a bath of anisotropic seed 
photons was studied both in the Thomson approximation and in the Klein- 
Nishina case and was used to produce optical and gamma-ray emission from 
the pulsars (see Morini 1981, 83 and reference therein). With reference to the 
hard X-ray properties of compact sources (AGNs and galactic black hole can- 
didates), Ghisellini et al. (1991) have derived (Thomson approximation) the 
angular distribution of the power and the spectrum of the hard IC radiation 
emitted by a relativistic electron located at the center of an hemispherical 
bowl, or cap, uniformly radiating soft X-rays (the emission region was as- 
sumed to be optically thin and the AIC scattering is taken to be isotropic 
in the electron rest frame). Dermer & Schlickeiser (1993) and Bottcher et al. 
(1997) have derived approximate solutions describing the AIC scattering of 
incoming photons from an accretion disc with the ultra-relativistic electrons 
in the relativistic jets of AGNs. 

Recently, Moskalenko & Strong (1999) have discussed the effect due to the 
anisotropic distribution of the incoming photons on the gamma-ray emission 
from interstellar photons IC scattered by cosmic-ray electrons in the Galaxy. 
This has been achieved by numerically integrating a redistribution function, 
in the Klein-Nishina regime, that was obtained by considering the scattered 
photons as unidirectionally emitted in the direction of the electrons. However, 
this approximation is strictly valid only in the ultra-relativistic case, which is 
the one of interest in the mentioned paper. 

The inverse Compton scattering of coUimated streaming particles with a sur- 
rounding bath of soft photons emitted by the disc-torus would produce high 
energy photons in AGNs and/or galactic X-ray binaries (Bednarek & Kirk 
1995, Protheroe et al. 1992). Furthermore the same mechanism has been used 
to model the gamma ray emission in the MeV blazar objects (Bednarek 1998). 
In order to better describe the emission pattern of all these astrophysical sit- 
uations general exact analytical expressions would be useful. 

The aim of the present paper is to investigate the AIC single scattering in the 
general case: from the trans-relativistic to the ultra-relativistic limit, both in 
the Klein-Nishina case and in the Thomson approximation. General formulae 
of astrophysical interest are derived. Specifically, we will consider the scatter- 
ing of a photon beam with an electron beam (Sect. 3) and with an isotropic 
distribution of electron momenta (Sect. 4). In Sect. 2 we define the basic ge- 
ometry of the scattering and obtain general relationships from the scattering 
between a monochromatic photon beam and an electron. In Sect. 4 we also 
compare our results with those of Bonometto et al.(1970) and Aharonian & 
Atonyan (1981) for the ultra-relativistic AIC and with the isotropic ultra- 
relativistic inverse Compton case (Blumenthal & Gould 1970) . 
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As an application, in Sect. 5 we use our general AIC equations of Sect. 4 to 
study the X-to-gamma ray spectrum of the powerful radio galaxies deriving 
from the AIC scattering of the nuclear IR-to-UV photons from a hidden 
quasar by the relativistic electrons in the radio lobes. In particular we argue 
that the detection of such emission with the new generation X-ray telescopes 
(Chandra, XMM) can provide important information about the magnetic field 
strength, the energy density, distribution and aging of the relativistic particles 
in the radio galaxies. 



2 IC scattering between a photon beam and a single electron 

We assume a scattering geometry such that the momenta of the incoming and 
scattered photon and of the electron are described by the four- vectors (Fig.l): 



P= -(1,0, 0,1) (1) 
c 

Pi-'jihk,) (2) 

Pe = mc7(l, Pcj) (3) 

e and ei being the incoming and scattered photon energies, kj, ej (j = 1, 2, 3) 

the unit vector components in spherical coordinates, m the electron mass and 
P — v/c the electron velocity. 

With these assumptions the incoming photon density is given by: 

n{e, Qph)de(fQph = 7^S{e - €o)S{cos9ph - l)de(PD.ph (4) 



eo being the assumed incoming photon energy, 5 the Dirac-function. 

It is well known that in the Compton scattering the tetra-vector components 
of the incident and scattered photons in the lab frame and in the electron 
frame (primed) are related by (Rybicki & Lightman 1979, Pozdnyakov et al. 
1983): 

ei=7(l + /9e;.fc;)e; (5) 
e' = 7(1 - Pes)e (6) 
e[ = e'(l + e'(l - cosesc)/mc^)-^ (7) 

where 7 is the Lorentz factor of the electron. 
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Fig. 1. The geometry of the scattering: the incoming monochromatic (energy e) 
photon beam propagates along the Z-axis and the photons are scattered to an energy 
ei in the direction {9sc^4'sc) by electrons with momenta in the {O^cj)) direction. 

In general the inverse Compton emissivity in the lab frame may be obtained 
by the manifold integral (Blumenthal & Gould 1970): 



j(^sc,ei) = j de d^Q!pf^d^Q,ed^ 



d^n'je', e, Vlph) d^n'scde[dt' 
de'd^nU d^Qscdeidt 



'ph 

d^a 



d'^Q'scde[ 



(8) 



where d^n', the photon number density in the electron frame, can be easily 
derived from the relativistic invariant dn/e (Blumenthal & Gould 1970). For 
a monochromatic (eo) photon beam moving along the Z-axis (Eq.4) it is: 



d^n'{e', fl'pfj^; e, Qph) _ n e' d'^Vtphde 
demQ^ ~ 2^1 d'^n'pi^de 



^fl^5{cosepn - l)5{e - 6o) (9) 



da is the differential Compton cross section, in our case the unpolarized Klein- 
Nishina cross section (Berestetskii et al. 1982): 
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1 + cos' 



\^e'{\-cosQ'sc)lmc^ 



(10) 



A^e(7) ^e) is the electron number per energy and solid angle. In this case it is 
formally written as: 

7V,(^, ^,)d^S^, = S^Qe - - lo)djd^^e (11) 



and cos9'gQ, the scattering angle in the electron frame, is derived by a Lorentz 
transformation of Eqs.(l-2) along the electron direction : 



The velocity of the electron that scatters a photon from energy e to ei is derived 
by combining Eqs.(5-7) and making use of the relativistic transformation of 
the angles {Pcjkj ^ 1): 



(3 = {{cjkj - ees/ei)^ + 

i-h 



e)TH(i--)(eA-— ) + 



6{(e,% - ^)2 + (- 

61 mc 



e 



(13) 



that in the Thomson approximation (i.e. je{l — k^) « mc^) becomes 

/3 ^ (14) 
eiCj/tj — £63 



The Compton emissivity is readily obtained by Eqs.(8-12) and (5~7) and well 
known relativistic transformations: 



j(^sc,ei) 



5(e-eo){(l + 



ei /C3 - 1 , 



r+[i + 



h-i . 

£1 ^3 - 1 



} (15) 
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where L = 1 — Poe^, Li = 1 — Poe'jkj, Pq is given by Eq.(13) with ej and 63 
replaced with the assumed electron coordinates (e° and 63 respectively) and 

ep, = -^{1 - ^(^3 - 1) - -L} (16) 
fJo 7omc^ ei 

In the Thomson approximation Eq.(15) becomes the much more simple: 
j(l^.c,e,)-^^5(e-eo){2 + 2-^ + (-^) } (17) 

with j3o given by Eq.(14) and 

e5A;,^l{l--L} (18) 



3 IC scattering between photon and electron beams 

As previously assumed the photons propagate along the Z-axis. In this Section 
the direction of the electron beam is defined by the coordinates 63 = Cg = 
cos 6^ and = 0°. 

The electron number density is formally given by: 

where and f{pf) are the number density and energy distribution of the 
electrons. 

In this case given eo, Ci, the electron direction Cj and the scattered photon 
direction kj, only one electron energy (7) satisfies the conservation laws. As 
a consequence the ^-function appearing in Eq.(15) can be written as: 

5{e - eo) = 5(7 - l)(3i''J (20) 

where 7 = (1 — /3^)~^/^ is given by Eq.(13) with eo instead of e, while the 
Jacobian of (3 with respect to e is: 

^ An alternative form of Eq.(17), but less useful for the purpose of this paper, has 
been recently obtained by Fargion et al. (1997). 
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Fig. 2. The geometry of the scattering as seen along the X-axis. The incoming 
photon beam and the scattering electron beam are shown together with the emission 
cones for 63 < eg and 63 > 63. 



J = {(e,% - ee,le,f + ((1 - k^)elmc^fY\2i5t-^{e,k, - ^) 

ei ei 

ejkj + 63-2^ 



-( 



} 



(21) 



2m(? $^ 

with $ = ((1 - k^)e/mc^Y + (e^^ - ee^/eiY - (1 - e/ei)^. 

The Compton emissivity is obtained by integrating Eq.(15) over the electron 
distribution (Eq.l9). By making use of Eqs. (20-21) and well known relativistic 
transformations, the integration yields: 



where f3o is given by Eq.(13) with the assumed electron coordinates e° and 
the incoming photon energy eg replacing the general ej and e respectively. J° 
is given by Eq.(21) with e'jkj, 63, and (3q replacing ejkj, 63 and (3, respectively, 
and where 
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^o = {l 



mc2 7o(l-^oe§) 



+{1 + 



(23) 



Eq.(22) can be used to calculate the emission due to AIC scattering between 
a streaming population of relativistic electrons and seed photons emitted by 
a small region (e.g., a star close to the jet or a hot spot in the accretion 
disk) in AGNs or galactic X-ray binaries. To describe with high accuracy 
much more common astrophysical situations the solution can also be easily 
obtained by numerically integrating over the distributions of the photon and 
electron momenta. 

In order to describe the AIC emitted power and spectrum and the changes 
introduced by the non ultra-relativistic scattering kinematics, we have to as- 
sume an energy distribution of the electron population. Let us assume a simple 
power law energy distribution, 7(7) = 7"'^, so that the results of this section 
can be compared with those in which the electron momenta are isotropically 
distributed (next section). It should be born in mind, however, that there are 
no known simple mechanisms providing a power law energy distribution in 
accelerated electron beams (Fermi mechanisms in general work with isotropic 
electrons momenta). 

The AIC emissivity is simply obtained by replacing 7(7 = 70) with 7^'' in 
Eq.(22). 

It can be noticed that the emitted power vanishes when 70 goes to infinity {S > 
1). Electrons of increasing energy are required to scatter the photons to a given 
energy ei at larger angles from the electron beam, up to a maximum angular 
distance for which (3 = 1. Therefore photons of energy ei are scattered in the 
line of sight only if the angular distance of the electron beam is sufficiently 
small. Since the scattering problem is symmetric with respect to the Z-axis, 
it is simpler to study the spatial distribution of the scattered radiation by 
assuming that the electron beam lies on the (Y-Z) plane. From Eq.(13), with 
sin(f)sc = sin(f)^^ = 1, and e = eo, one finds that the limiting condition (3o — ^ 
is obtained at two extremes in (Fig. 2): 



M±) = e°[(l-6o(l-e°)/ei]±{(l-(e°)2) 

•(1 - (1 - ^(1 - emy/' (24) 



Similarly, by releasing the condition sin (psc = 1, one obtains the extremes 
in (psc by imposing the condition k3{+) = /c3(— ). If 0^ is large enough, the 
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Fig. 3. The emitted power as a function of the scattering angle for the IC scattering 
between a photon and an electron beam in the Thomson approximation (solid lines) 
and with the Klein-Nishina equations (dotted lines). The curves are normalized to 
the maximum in the Thomson approximation. The electron beam is positioned 
at 0e(i'ad)=2.1O and = 7i"/2. The upper and lower (solid and dotted) curves 
correspond respectively to 6= 2.5 and 5.0, while the assumed initial and scattered 
photon energies are eo = lO^eV and ei = lOMeV. 

emission cone may comprise the Z-axis. The angle of the electron velocity at 
which this happens is derived from Eq.(24) by setting /c3(— ) = —1: 



For 63 > 63 the scattered directions on the (Y-Z) plane are comprised in 
the interval A;3(+) < < A;3(— ) and (f)sc — 7r/2, while for 63 < 63 the 
scattered directions are bounded by — 1 < A;3 < A;3(+) and (psc = 7r/2 and 
-1< A;3 < ks{-) and (f>sc = 37r/2 (Fig.2). 

In the Thomson limit Eq.(22) becomes: 




-l + ep/ei 
1 + eo/ei 



(25) 




(1 



7o 



ih - 1) 



1 



2 + 



10 




Log e,(eV) 



Fig. 4. The calculated spectrum from IC scattering between a photon and an 
electron beam is shown in arbitrary units for different scattering angles both in the 
Thomson approximation (solid line) and in the Klein-Nishina scattering (dotted 
line). We have assumed 5=2.5 and an incoming photon energy eo = 10^ cV. The 
assumed scattering angles from the bottom to the top of the diagram are: 6sc 
(rad)= 2.20, 2.15, 2.13, 2.11 with the electron beam positioned at 0g(rad)=2.1O and 
0O = 7r/2. 

^h^^ eoelil^-eiri^^'"^ ^^^^ 



where 70 and Pq are given by Eq.(14) with ej replaced by the assumed electron 
beam coordinates e° and with e replaced by eo- Furthermore it can be readily 
shown that in the Thomson approximation the emitted pattern has the same 
boundaries of the Klein-Nishina case. 

An example of the emitted power as a function of the scattering angle, both 
in the Thomson approximation and in the Klein-Nishina case, is represented 
in Fig. 3 for two different values of the electron spectral index. Relatively more 
power is channeled at larger scattering angles with decreasing S. The broad- 
ening in the emitted power distribution in the case of low values of 6 is due to 
the contribution of the high energy electrons that emit along relatively large 
values of 
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Due to the Compton recoil the energies of the electrons that are required to 
scatter incoming photons (eg) up to an energy ei in the Klein-Nishina regime 
are larger than those in the Thomson approximation. As a consequence the 
emitted power in the Klein-Nishina regime is more depressed in the case of a 
steeper energy distribution of the electrons (as it is 5 = 5 in Fig. 3). 

The emitted power per unit frequency and solid angle as a function of the 
emitted energy is represented in Fig. 4 for different scattering angles both in the 
Thomson approximation and in the Klein-Nishina case. The very sharp cut- 
off in the emitted spectrum is simply due to the kinematics of the scattering 
and only weakly depends on the electron differential spectrum. Such feature 
represents the portion of the emitted spectrum contributed by the electrons 
with (3 approaching to unity and with the emitted photons moving along the 
boundaries of the scattered cones. 

In Fig. 4 it is seen that the difference between the Thomson approximation 
and the Klein-Nishina scattering becomes larger at small scattering angles 
(with respect to the direction of the electron beam) and increases with the 
energy of the scattered photons. Actually, for each scattering configuration this 
difference follows a quite complicated trend being smaller both at low and at 
high emitted energies and larger at intermediate energies. The trend is due to 
the difference between the Thomson and the KN cross sections but, because 
of the Compton recoil, it is also sensitive to the electron energy distribution 
(see Brunetti 1998 for further details). 

A relevant example is the Thomson spectrum in the direction of the electron 
beam (e^/cj = 1 and = Cg). In this case from Eq.(26) the Compton emissivity 
is: 



^ ei — ep ^2 2eieo — (1 + e3)eo ^2-j 
^1 — 6360 6160 

,0\,2 



{2eieo-(l + e^)e^}— (27) 



In the ultra-relativistic limit Eq.(27) can be expanded in series for ei >> eo 
to get the much more simple: 

j(e°,60 - crlnK,2'^{l - e°)^(^)-^ (28) 



Now, as one can notice from Eqs.(5-7) if 9sc << 1 then 70 >> 1, however, 
the ultra-relativistic condition ei >> eo is not satisfied. As a consequence 
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Fig. 5. The calculated spectral index (j oc of the spectrum from a scattering 

between a monochromatic photon beam and an electron beam with a power law 
energy distribution {S = 2.5). In the scattering the photons are emitted in the 
direction of the electron beam. Prom the bottom of the diagram wc have reported 
the spectral index by assuming respectively ei/eo=50, 100, 200, 1000. The dotted 
line is the spectral index in the ultra-relativistic case (Eq.28). 

although in the case 6° ~ 1 is 70 — > 1, the spectrum differs from the so-called 
ultra-relativistic case as it can be noticed from Eqs. (27-28). 

The spectrum per unit solid angle in the mildly-relativistic case (Eq.27) is 
rather comphcated. In Fig. 5 its slope is shown as a function of the electron 
beam direction and of ei/eo; it is also compared with the ultra-relativistic 
value (Eq.28). 

The power emitted by the electron beam over Air can be compared with the 
formulae giving the emitted power in the case of isotropic inverse Compton 
scattering given in the hterature (Jones 1968, Blumenthal & Gould 1970, 
Rybicki & Lightman 1979). In the ultra-relativistic case the integration of 
Eq.(26) over flsc yields: 




o ,61, <5-l S+3 , n. c^+l 

-lnK,{^)- — 2—{\-e%) — 



(5^ + 4(5+11 



(29) 



(5 + l)(5 + 3)(5 + 5) 
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Log (e,/e) 

Fig. 6. The calculated spectral index of the spectrum integrated over 47r from a 
scattering between a monochromatic photon beam and an electron beam with a 
power law energy distribution {5 = 2.5) is shown as a function of the ratio ei/eo- 
From the top to the bottom of the diagram the electron beam direction is 6^= 0.5, 
1.5, 2.0, 2.5, 3.0 (rad), while is taken (pg = Tr/2. The spectral index is computed 
with 5 eV step (eo=l eV); in the ultra-relativistic limit it should be a = 0.75. 

the spectrum is a power law as in Baylis et al.(1967). In the mildly-relativistic 
case the spectrum is harder for large angles between the incoming photon and 
the electron beam, while it is softer at small angles (Fig.6). 



4 IC scattering between a photon beam and an isotropic electron 
distribution 

In the case of isotropic distribution the electron number density per energy 
and solid angle is: 

7V,(7, Q,)d^d''ne = {^n)-'KJ{^)d^d^n, (30) 

The emissivity is given by integrating Eq.(15) over the angular and energy 
distribution of the scattering electrons. The integral over can be performed 
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by following the method given in Aharonian & Atoyan (1981); one finds: 




, (7 — ei/mc^)ei/m(? + ^e^lrxK? + exe^k-^ j {rn(?Y 
^ [(7 - ei/mc2)2(l - k^f + 1 - A;i]3/2 

(7 + tQ/'mc^)eQ/'mc^ + ^ei/mc^ - eie^k^ / (mc^)'^ 
[(7 + eo/mc2)2(l - ^3)2 + 1 - A;i]3/2 ^ 

2~(mc^)^}n(i7 



6061 



(31) 



where the function 



7^(eo,el,7,A;3) = [(7 - ei/mc=^)'(l - A;3)' + 1 - fc^] 
-[(^ + eo/mc'f{l-ksf + l-kl] 




(32) 



The lower limit of the integral Eq.(31) is given by the minimum energy of 
the electrons that can scatter an incident photon (eo) to an energy ei in the 
direction ^sc- 

The minimum velocity of the electrons is obtained by minimizing Eq.(13) 
with respect to the electron directions {cj, i.e. 9e,4>e)- Given that the distri- 
bution of the scattering electrons is isotropic, it is equivalent to describe the 
emitted spectrum with reference to any plane containing the incoming pho- 
ton momenta (the Z-axis). For simplicity of notation we study the scattering 
properties on the (Y-Z) plane, so that ki = 0, /c2 = (1 — ^3)^^^, and the elec- 
trons having the lowest energy to make the scattering are those moving on 
the (Y-Z) plane i.e. with 0e — '^/'^ (i-e. ei — 0). 

We find that Eq.(13) is minimized by electrons having directions : 



with e3(-|-) for > eo/ei and e3(— ) for k^ < eo/ei. By combining Eqs.(13) 
and (33) one finds Pmin and 'jmin = {1 ~ Pmin}^^^"^- should be noticed 
that Eq.(33) holds whatever the electron and the photon energy may be, i.e. 
both in the Thomson approximation and in the Klein-Nishina case, from the 
ultra-relativistic limit down to the trans-relativistic case. 



e3(±) = ±{1 + 



{ks - eo/ei) 




(33) 
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In the Thomson approximation, from Eqs.(14) and (33), we find: 



that becomes the simplest 7mm ^ (ei/2eo(l — ^3))^''^ in the uhra-relativistic 



case (Bonometto et ah 1970), and the weU known ^min = l/Syei/eo in the 
ultra-relativistic isotropic Thomson scattering (taking = —1). 



4.1 The Thomson approximation 



In the literature it is usually assumed a momentum-power law energy distri- 
bution of the relativistic electrons /{j) = since it is that expected 
from Fermi acceleration mechanisms. It deviates from a simple energy-power 
law in the non ultra-relativistic case. 

In this section our aim is to derive AIC equations from the trans-relativistic 
to the ultra-relativistic regime and to compare our findings with the ultra- 
relativistic equations in the literature. As a consequence, in order to evidence 
the changes in the AIC emitted spectrum and pattern introduced only by the 
non ultra-relativistic scattering kinematics, in this section, we adopt a simple 
energy-power law distribution /(7) = O'^^ down to trans-relativistic energies. 
The case of a general 7(7) is discussed in Appendix A where a simple semi- 
analytical solution is given. The case of the distribution f{j) = (3~^{(5j)~^ is 
discussed in Appendix B and is compared with the results of this section. 

From Eq.(31) the integration in the Thomson case yields: 



n{ 



2X0 



4 ^60^ ^((a)2_2A;3^i + l). 



2(1 - A;3)-^(l + ^)J3/2 + 

[1 - M-[(l + ^)(3A;3 -h- |(^ + (^)^)]X5/2 
5 



'''' ' ")(3A:3-p-^(^ + (^)^^ 
+ o(l-M-'[3(l + -) + (- + (-m/2} 



(35) 



where 

T = — V + i) ■ ,-{^+2»+l} /of^N 

° v^^r^((5 + 2i + l)r(i + l)^"^'" ^ ^ 
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Fig. 7. The normalized emitted power as a function of the scattering angle for dif- 
ferent values of 5: 2.0 (dashed line), 2.5 (dotted line), 4.0 (solid line). In the compu- 
tation it has been assumed eo=0.1 eV and ei=l keV i.e. Thomson ultra-relativistic 
case. Relatively more power is emitted al large 9sc with increasing 6. 
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I min ks 



gi,mT{m + l)T{t + 1) 



(3J 



I- 



7/2 



^ (m3 + X + fm + f)r(m+i)r(_^ 



where 6*^3 = (/cg - 1)/(1 - ksY, Qi^rn = 6 + 2{i + m) + 5, and 7™n is given by 
Eq.(34).Q 



^ It would also be noticed that the ratios r(z[m] -|- l/2)/r(i[m] -|- 1) in Eqs. (37-39) 
are also equivalent to 4'^^'^^B{i[m] + l/2,i[m] + l/2)/7r^/'^; B being the Beta function. 
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A solution in a simple closed form for Eqs. (35-39) that holds down to the 
mildly-relativistic case (for ei/eo > 15 — 20 with < 5% approximation) is: 

j(^3,6i) Ky,cn{l - h)"^ ' {A (5) + fMS,h)} (40) 



where 



and 



A 5^ + 4(5+11 

^'(*) = (^+l)(^ + 3)(^ + 5) 



, , , , , (6^ + 86^ + U6'^ + 1285 + 139) A;3 - 35^ + 5^ + 355 + 31 



This very simple formula could be very useful in astrophysical problems where 
non ultra-relativistic electrons are involved in the scattering process. 

By considering Eqs. (36-39) in the ultra-relativistic limit ei >> eg, after simple 
algebraic manipulation from Eq.(35) (or simply from Eqs. (40-42)) we find the 
emissivity in the ultra-relativistic case: 

i(k. e.) - K ,2,^il^MZ^(^!±i^±Il) fiLr^ f43) 



that is equivalent to the result given in Bonometto et al.(1970). Furthermore 
by integrating Eq.(43) over Qe (it should be noticed that due to the isotropic 
distribution of the electron momenta it is equivalent to integrate over flsc) 
one finds the well known solution in the isotropic ultra-relativistic Thomson 
scattering case (Blumenthal & Gould 1970). 

The emitted power is obtained by Eq.(35) in the general case or by the much 
more simple Eq.(43) in the ultra-relativistic case; in the ultra-relativistic case, 
it is shown in Fig. 7. As expected the scattered power has a maximum at 
dsc = 7^ and goes to zero for small scattering angles. The electrons that are 
mainly responsible for the scattered power at a given ei are those with energies 
close to the minimum energy that increases with the scattering angle. It follows 
that by increasing S there are relatively fewer electrons at higher energies and, 
consequently, the emitted power goes to zero more rapidly toward smaller 
scattering angles. 

In order to compare our general results with the AIC formulae in the literature 
we give in Fig. (8) the ratios between the AIC emission (Eq.35) and that of the 
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Fig. 8. The ratios between the general AIC and the ultra-relativistic (j^^*^') emitted 
power as a function of the scattering angle {9sc) are given for €i/e= 6, 16, 100. In 
the calculation it is assumed 7(7) = 7'"'^ with 5= 2 (solid lines), 2.5 (dotted lines), 
3 (small-dashed lines), and 3.5 (large-dashed lines). All the calculated patterns are 
normalized at 9sc = 180". One finds that the percentage of power emitted at small 
Osc increases in the mildly and trans-relativistic case. 

ultra-relativistic case (Eq.43) as a function of the scattering angle for different 
electron energy distributions and energies; the differences are > 30% also in 
the mildly-relativistic case (i.e. ei/e ~ 20 — 30). The main effect due to the 
non-relativistic kinematics is that the percentage of radiation emitted at large 
scattering angles decreases with decreasing ei/e. 

The emitted spectrum for different scattering angles and 5=2.5 is shown 
in Fig.9. As already noticed (Eqs.35 and 42), at variance with the ultra- 
relativistic standard equations in the literature, the AIC spectral index cal- 
culated with our general equations is a function of the scattering angle and 
of the energy. With decreasing ei/e the spectrum flattens at large 9sc and 
steepens at small 9sc\ this causes the progressive decrease of the percentage 
of radiation emitted at large 6sc (Fig-8). The difference from the standard 
ultra-relativistic spectrum (that is a — 0.75), vanishes at higher energies due 
to the relativistic aberration. 
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Fig. 9. The calculated spectral index (j oc z/~") computed with a 5 eV step with 
eo=0.1 eV) in the Thomson regime as a function of ei/eo for several values of the 
scattering angle. Starting from the top the scattering angles are: 05C'(rad)= 0.8, 1.2, 
1.5, 2.8, 3.1. We have assumed S=2.b and eo=0.1 eV. At lower energies and large 
scattering angles the spectrum is flatter than in the highly relativistic isotropic case 
(a = 0.75), while it is steeper at small scattering angles. The differences become 
vanishing with increasing the emitting frequency. 



4-2 The Klein-Nishina case 



It is well known that the emitted IC spectrum, in the case of an isotropic 
distribution of the seed photons momenta, typically steepens at an energy 
{^kn) such that e^^y ~ mc^. 

In the case of the AIC scattering ekn is expected to depend on the scattered 
directions. Eq.(31) can be used to investigate the Klein-Nishina steepening 
at different scattering angles both in the ultra-relativistic and in the general 
case. This point is illustrated in Fig. 10 where we have plotted the AIC spectra 
for different values of 6sc by assuming a power law energy distribution of 
the scattering electrons (/(7) = 7^^, 5=2.5 and eo=l kcV). The behaviour 
in Fig. 10 is readily understood by keeping in mind the kinematics of the 
Compton scattering (Eqs. (5) and (7)). One has that ckn ~ ~ jmc^ 

and that 7, the typical (or minimum) energy of the electron required to scatter 
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Fig. 10. The calculated spectral index (computed with 5 MeV step; j oc in the 
Klein-Nishina regime for an isotropic distribution of the relativistic electrons as a 
function of the emitted energy for different values of the scattering angles. We have 
assumed 5=2.5 and a monochromatic incident photon beam (eo = 10^ eV). Starting 
from the bottom of the diagram the STicccssive values of the scattering angles are: 
9sc{^SLd)= 0.2, 0.5, 1.0, 1.5, 2.0, 2.8. The Klcin-Nishina spectrum steepens at higher 
energies with decreasing scattering angle. 

a photon (eo) to ei increases for large 9sc- As a consequence the KN-phase 
starts at smaller emitted energies with increasing 9sc- On the other hand, we 
have shown that in the Thomson approximation the emitted power is greatly 
enhanced at large 6sc (Fig. 7). These two effects can largely compensate each 
other. This is clearly illustrated in Fig. 11 where we have plotted the emitted 
power as a function of the scattering angle (with 6=2.5 and eo=l keV) for 
several values of the scattered photon energy ei. When the KN regime is 
well established (larger ei) the scattered power is much more isotropically 
distributed (a sort of KN isotropization) . 

As shown by Aharonian & Atoyan (1981), the shape of the AIC emitted spec- 
trum in the ultra-relativistic case is identical to the isotropic case. Aharonian 
& Atoyan (1981) Equation can be obtained by integrating Eq.(31) over the 
electron spectrum in the ultra-relativistic limit. 

However, since the KN regime may become important at relatively low energies 
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Fig. 11. The normalized emitted power as a function of the scattering angle in the 
Klein-Nishina regime is plotted for three different values of the emitted energy: 
ei= 10^ eV (solid line), 10^ eV (dotted line), 10^^ eV (dashed line). We have as- 
sumed 5 = 2.5, and an incoming photon energy eo= 10^ eV. The emitted power 
is distributed much more isotropically (KN-phase isotropization) at larger emitted 
energies. 

of the scattering electrons, both the non ultra-relativistic effects previously 
discussed (Fig.9) and the classical KN steepening may show up in the emitted 
spectrum. This happens to be in the energy range that simultaneously satisfies 
eieo(l — /cs) > 10~^(mc^)^ (for which the spectrum becomes sensitive to the 
KN effects) and ei/eo < 10^'^ (i.e. non ultra-relativistic). It is illustrated in 
Fig. 12 where we have plotted the spectral index as a function of the scattered 
photon energy for two representative values ofOsc (with S=2.5 and eo=l keV). 
Our calculations show that at low energies the spectral index is much steep in 
the case of the smaller 6sc and starts to converge toward the classical ultra- 
relativistic limit (a=0.75 in this example) but the progressive KN steepening 
of the spectra takes over before the ultra-relativistic regime is estabhshed. 
When the KN regime is well established the spectral index is much steep in 
the case of the large 9sc- 
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Fig. 12. The calculated spectral index (j oc v~°^ computed with a 10^ cV step) in 
the Klein-Nishina regime as a function of the emitted frequency for two different 
values of the scattering angle: Bsc= 1-2 (dotted line), 2.4 (dot-dashed line). We have 
assumed 5=2.5 and an incident photon energy eo =10^ eV. In the diagram both the 
departure from the classical isotropic spectral index at low photon energies and the 
spectral steepening due to the Klein-Nishina effect are shown. 

5 Application to the FR II radio galeixies 



By adopting the unification scheme linking radio loud quasars and FR 11 
radio galaxies (Barthel 1989), it has been shown that the AlC scattering of 
the nuclear photons by the relativistic electrons in the lobes may produce 
detectable X-ray fluxes well in excess of those calculated by the IC of the 
CMB photons (Brunetti et al.l997). 

Roughly speaking, the inverse Compton emissivity is oc K^LqsoI ^ R and 
Lqso being the distance from the hidden quasar and the far IR-UV quasar 
luminosity respectively. By constraining Lqso, a comparison of the X-ray 
properties with the radio synchrotron emission can provide important infor- 
mation on the magnetic field strength and relativistic particle distribution in 
the lobes. 

In a recent ROSAT HRI observation of the FR II radio galaxy 3C 219 Brunetti 
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et al.(1999) have found extended non-thermal emission from the lobes most 
likely due to AIC of nuclear photons by relativistic electrons. It was found 
that the average magnetic field strength of 3C 219 is a factor ~3 lower than 
that computed under equipartition assumption. 

In order to illustrate the possibility of obtaining detailed physical information 
from X-ray observations, which will be possible with the future Chandra and 
XMM missions, in this Section we apply the general AIC equations of Sect. 4 
to a specific model. 

5. 1 Calculation of the emissivity 

We assume that the far IR-optical spectrum of a typical quasar hidden in 
a powerful radio galaxy is represented by four power laws (-^,7") with a — 
0.2, 0.9, 1.7, and 0.6 respectively in the intervals 100-50, 50-6, 6-0.65, and 
0.65-0.35 (as in Brunetti et al.l997). Since the UV photons may be important 
in the IC production of the hard X and gamma-rays, in addition we assume 
a UV spectrum simply modelled by a single power law with a = 0.0 in the 
0.35-0.03/im band, consistently with Walter & Fink (1994) findings. 

We further assume that the relativistic electrons are injected in the radio 
volume, approximated with a prolate ellipsoid of 100 x 50 kpc semiaxis, with 
an energy differential spectrum Ne{j) = Kej^^ (the calculated AIC fluxes 
are typically contributed by electrons by 7 > 5 — 10) and that the magnetic 
fields, of average constant intensity B, and particles momenta are randomly 
distributed on a sufficiently small scale. 

The time evolution of the particle spectrum, injected at the hot spots and no 
reacceleration, due to synchrotron and Compton losses may be described by 
(Jaffe & Perola 1973): 

N,{-f) oc 7-'5(l - ^)^-2; for -f<-fb (44) 

76 

being A^e(7) = for 7 > 76 and % the break energy (Kardashev 1962). This 
is crucial in order to compare the predictions with detailed spatially resolved 
X-ray spectroscopy. Deep radio observations of the hot spot regions of the 
radio galaxy Cyg A have suggested the possibility of a low energy cut off, or a 
turn over in the particle distribution (Carilli et al. 1991). We take into account 
this possibility by applying an exponential cut off to the power law spectrum 
below 7c. 

In general, is expected to be a function of the position R in the radio galaxy. 
In a very simple model, in which the radio hot spots separate with constant 
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Fig. 13. The calculated break energy (in units of mc^) is reported as a function of 
the distance from the nucleus Rkpc for different magnetic field strengths i?=l-10~^ 
(A) and 5- 10"^ G (B). A distance of the hot spot from the nucleus £)=100 kpc, 
a velocity of separation 77=0.04 and a rcdshift z=0.15 have been assumed. The 
luminosity of the hidden quasar from the bottom to the top of the diagram (solid 
lines) are: 5 and 1 • 10^^ erg s~^. The dotted lines represent the prediction without 
the hidden quasar. 

velocity (rjc) and the particles remain approximatively in the same place in 
which they are injected, the age of the relativistic particles ist — (D — R)/ {rjc), 
D being the distance of the hot spot (for which t—0) from the nucleus. 

At variance with the previous literature, in addition to the standard syn- 
chrotron and Compton (with the CMB photons) losses, we consider also the 
radiative losses due to the AIC scattering of the nuclear photons. 

With these assumptions the break energy of the time/spatially evolved elec- 
tron population in the radio galaxy is: 

7.(i?..c) = 5-10-^^ ^2^2 + ^2^(1 + + 7 . (45) 

J^kpc J^kpc <J ^kpc 

where Bic is the equivalent magnetic field associated to the CMB and Lq^q is 
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the luminosity of the hidden quasar in units of lO^^erg s~^. The AIC scatter- 
ing with the nuclear photons dominates the radiative losses at relative small 
distance from the nucleus. This distance depends on the quasar luminosity 
and on the magnetic field strength in the lobes (Fig. 13). 

Due to the decrease of the break energy with decreasing distance from the 
nucleus, the innermost regions of the radio volume contribute less to the ra- 
diation spectrum of the radio galaxies at high frequencies. Spatially resolved 
spectroscopy in the X-ray domain can provide important information about 
7(,(i?), especially in the innermost regions of the radio galaxies. The break en- 
ergy can also be constrained by deep radio spectroscopy. However, it should be 
noticed that in the innermost parts, where 7fe(-R) < 10^, the expected repen- 
tine decrease of the radio brightness makes these investigations very difficult. 
If the source is well resolved both in the radio and in the X-ray bands, then 
Lqso-i B and rj can in principle be estimated by comparing the radio and the 
X-ray spectral indices with the model predictions. 

In Fig. 14, for a given radio galaxy model, we have reported the spectral index 
from the AIC scattering of the nuclear photons at different distances from the 
nucleus. Fig. 14 is obtained by integrating Eq.(31) over the electron distribu- 
tion (Eq.44), with 7^ given by Eq.(45), and over the assumed spectrum of the 
nuclear photons. For simplicity, the axis of the quasar radiation cone, of half 
opening angle 45°, is assumed to be coincident with that of the radio source 
and to lie on the plane of the sky. In order to give theoretical predictions to be 
compared with the future X-ray observations, a typical radio galaxy of total 
size 200 kpc is placed at a distance (z=0.15) such that its angular size is much 
larger than the Chandra PSF. Since the 0.2-10 keV luminosity of a strong 
FR II radio galaxy is expected to be close to ~ 3 • 10^'^erg s~^, the estimated 
Chandra count rate from the two regions indicated in Fig. 14 is ~ 10~^ counts 
s-i {Ho=7b km s"^ Mpc'S go=0.0). 

The complex spectral index behaviour is due to the interplay of the IR and 
UV bumps with the time/spatially evolved particle energy distribution. The 
X-ray spectral index of the innermost regions of the radio galaxies is signifi- 
cantly steeper than those pertaining to the extended regions and to the radio 
synchrotron spectral index of the high brightness lobes. 

We have also computed the X-ray spectral index and emitted power by im- 
posing a low energy cut-off in the particle energy distribution (dashed lines) . 
In this case the soft X-ray spectral index (0.1 - 1 keV) is expected to be 
considerably harder than the synchrotron spectral index (=0.75). As a con- 
sequence detailed spatially resolved spectroscopy provided by Chandra may 
yield unprecedented information about the energy distribution of the relativis- 
tic particles in the low energy portion of the spectrum. 
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Fig. 14. The calculated spectral index (j oc due to the AIC scattering of incom- 
ing IR-UV photons from a typical quasar by a time-cvolvcd isotropic population 
of relativistic electrons N{'y) oc with no energy cut-off (solid line) and with 

a low energy cut-off at 7 = 20 (dashed line). The assumed parameters are those 
of Fig. 13 (case A) with the luminosity of the hidden quasar 5 • lO^'^crg s-^. The 
radio volume is approximated with a prolate ellipsoid of 100 x 50 kpc semiaxis. The 
spectra are calculated for two regions comprised between 2-7 arcsec (2) and 20-25 
arcsec (1) from the nucleus ( 1 arcsec = 2.36 kpc; Ho=75 km s~^ Mpc~^, qo=0.0). 
The synchrotron spectral index of the high brightness radio lobes is 0.75. 



If the radio galaxy is inclined with respect to the plane of the sky the nuclear 
photons in the far lobe are on average scattered toward the observer at angles 
larger than those in the near lobe. Since in the case of AIC scattering the 
typical 7 of the scattering electrons depends on the scattering angle (Eq.34), 
one finds that in the far lobe the nuclear photons are scattered at a given ei by 
less energetic electrons than in the near lobe. As a consequence, depending on 
the inclination of the radio galaxy on the sky plane, all the spectral features 
in Fig. 14 would appear shifted at lower energies in the case of the near lobe 
and at larger energies in the case of the far lobe. Typically this results in a 
spectrum of the far lobe significantly harder than that of the near lobe. 



In the case of distant (or faint) radio galaxies spatially resolved spectroscopy 
will not be possible. However, due to their large effective area and spectral 
resolution, Chandra and XMM could be used to study the integrated spec- 
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Fig. 15. The calculated spectral index (j oc due to the AIC scattering of incom- 
ing IR-UV photons from a typical quasar by an unresolved source. The assumed 
magnetic field strength are: 5 (solid lines) and 30-10"^ G (dashed lines). It is shown 
for different hidden quasar luminosities, from the top to the bottom of the diagram 
one has: 10 and 2 -10^^ erg s~^. A distance of the hot spot from the nucleus £>=100 
kpc, a velocity of separation 77=0.02 and a redshift 2;=0.15 have been assumed. The 
radio volume is approximated with a 100 x 50 kpc prolate ellipsoid. The injection 
index is 5=2.5 so that the radio spectral index of the high brightness radio lobes is 
a=0.75. 

trum. 

The total emitted power is obtained by integrating the emissivity over the 
emitting volume. The contribution from the IC scattering of the CMB photons 
should be taken into account in the case of high redshift and/or very extended 
radio galaxies. The radio volume can be modelled in detail when one deals 
with a specific radio galaxy. In Figs. 14 and 15 the calculated AIC spectral 
index and emitted power (in arbitrary units) from powerful radio galaxies are 
shown from the soft X-rays to the gamma-rays for a number of parameters. 
The spectral index from unresolved radio galaxy is close to the radio spectral 
index in the soft X-ray band, being steeper at higher energies. In the case of a 
powerful radio galaxy at a redshift 2;=0.5 with radio spectral index a= 0.75 the 
expected spectral index in the XMM band would be (Xx= 0.75-0.90 depending 
on the parameters (i.e. quasar luminosity, age, magnetic field strength). If the 
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Fig. 16. The calculated AIC emitted power per unit solid angle (normalized at 100 
eV and in arbitrary units) due to the scattering of incoming IR-UV photons from 
a typical quasar by an unresolved radio galaxy is shown for different luminosities of 

the hidden quasar: 2 (top curves) and 10 -10^^ erg s~^ (bottom curves). The assumed 
magnetic fields strength are 5 (solid lines) and 30-10"^ G (dashed lines). The other 
parameters are the same of Fig. 15. 

radio galaxy is inclined with respect to plane of the sky the spectrum from the 
far lobe would be harder than that of the near lobe. Since the far lobe is also 
expected to be the most luminous (Brunetti et al.l997), then the integrated 
spectrum from an unresolved radio galaxy would be harder with increasing 
inclination. 

Prom Fig. 16 we also notice that a large fraction of the energy emitted by the 
AIC scattering of the far IR-UV nuclear photons is channeled in the gamma- 
ray band. This will provide a new class of extragalactic gamma-ray sources 
to be revealed by future high energy experiments such as GLAST. 



6 Conclusions 

The anisotropic inverse Compton (AIC) scattering problem has been solved 
from the trans-relativistic to the ultra-relativistic regime, without introducing 
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any approximation. 



The equations are integrated over the energies of the scattering electrons to 
give the spectrum from the AIC scattering of an unpolarized monochromatic 
photon beam by a beam of unpolarized relativistic electrons. In the case of 
a power law energy distribution of the electrons, we find that the spectrum 
is inverted peaking toward high energies before being sharply cut off at the 
emitted energies at which the highest energetic electrons (/3 — > 1) are involved 
in the scattering process. The emitted power and spectrum in the Thomson 
approximation is compared with those in the Klein-Nishina regime. In the 
ultra-rclativistic Thomson approximation the radiation scattered in the direc- 
tion of the electron beam has a power law spectrum with slope a = (5 — 3)/2. 
We find that in the mildly-relativistic case it is much harder at large scattering 
angles and much softer at small scattering angles. 

The AIC scattering between unpolarized isotropic electrons and a monochro- 
matic photon beam is discussed from the trans-relativistic to the ultra-relativistic 
case. In agreement with the previous studies, we find that, toward the ultra- 
relativistic limit, the distribution of the scattered power is highly anisotropic 
for the astrophysically common case (where the differential electron spectral 
index S > 2) and that the spectrum has the same shape of the isotropic ultra- 
relativistic case. In addition to previous studies we find that in the mildly- 
relativistic case the spectrum emitted at large scattering angles is considerably 
harder than the classical ultra-relativistic one being softer at small scattering 
angles. Furthermore, we show that the emitted power in the Klein-Nishina 
phase is distributed much more isotropically with increasing energy of the 
scattered photons. We have also derived simple formulae that give the emis- 
sivity from the ultra-relativistic case down to ei/e ~ 10. 

In Sect. 5 we have applied the general AIC equations derived in Sect. 4 to the 
calculation of the emitted power and spectrum from the scattering of the far 
IR-UV nuclear photons from a hidden quasar by the relativistic electrons in 
the lobes of powerful FR II radio galaxies. In our simplified model we have 
taken into account the effect of the radiation losses on the energy distribution 
of the relativistic electrons. Due to the losses the calculated spectrum steepens 
with decreasing distance from the nucleus where the AIC brightness increases. 
This may provide an important test of the model to be performed by future 
X-ray satellites such as Chandra and, conversely, it may provide extremely 
important information on the magnetic field strengths and on the time evo- 
lution of the relativistic particles in the FR II radio galaxies. In the case of 
resolved radio galaxies inclined with respect to the plane of the sky, due to the 
anisotropic illumination from the nuclear photons, we predict that the X-ray 
spectrum of the far lobe is on the average harder than that of the near one; 
since the far lobe is expected to be the more luminous, also the spectrum of 
unresolved sources will be harder with increasing inclination. 
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7 Appendix A: Thomson AIC scattering by isotropic electrons 
with a general energy distribution function 



In Section 4.1 we have derived the Thomson AIC emissivity by an isotropic 
electron population with energy distribution /(7) = 7"*^ (Eqs. 35-39). In sev- 
eral astrophysical situations f{'f) may be more complicated than a simple 
infinite power law (e.g. Sect. 5). In these cases the Thomson AIC emissivity is 
given by: 



j(/c3,ei) ~ — r^(-)M- 



4 ^60^ ^((a)2_2A:3a + i), 



2(l-fc3)2(l + -H)X3/2 + 
(1- ^3)^(1 + 7)(3^3-|)-|(5i + (^)^)]%2 

+|(1 - hn3{l + ^) + (^ + C-^my2} (46) 
z ei eo ei 



with 



Imin 

= I (^^^!^{7^(i - hf + 1 - kir'^'d^ (48) 

^5/2 = / '1^2^)1/2 (7^(1 - ^3)^ + 1 - klr^^'d^ (49) 

%2 = / (^^^^{7^(1 - + 1 - klr^'d^ (50) 
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Fig. 17. The ratio between the angular distribution of the scattered radiation from 
Eqs.(35 and 51-54) and from Eqs. (35-39) (i.e. simple power law in 7; j*) is shown 
as a function of the scattering angle. In the calculation 6=2.5 has been assumed; the 
distributions are given for ei/e=16 (solid line), 10 (dotted line), 5 (small-dashed 
line), and 4 (large-dashed line). Due to the decreasing energy of the scattering 
electrons, the differences between the two calculations are larger for large scattering 
angles. 

which, in general, can be readily numerically calculated with 7^j„ given by 
Eq.(34). 



8 Appendix B: Thomson AIC scattering by isotropic electrons 
with /(p) = 



It is well known that Fermi acceleration mechanisms lead to an energy dis- 
tribution of the relativistic particles that is a power law in momentum and 
not in 7. As a consequence, in some astrophysical situations in which trans 
and mildly-relativistic AIC scattering produces radiation at relatively large 
scattering angles (typically > 50 — 60") one should use /(7) = /3~^ {/S'y)'^ . 
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In this case the emissivity is calculated from Eqs. (46-50) with /(7) = P ^(/^t) 
In analitycal form it is given by Eq.(35) with Eqs. (36-39) replaced with: 



r = 1 V r(i + 1 + 5/2) -{,+2^+1} 

° r(V2 + i)^(5 + 2i + i)r(i + i)^"^^" ^ ^ 



^ _ 2 - (m+l)r(m + l)r(i + l + V2) -^,^-2} .... 



J 4 - (m^ + 2m + f )r(m + |)r(i + 1 + ^2) ,„ 

3V^,£lo ^,^r(m + i)r(i + i)r(5/2 + 1) ^ ^ 



x,/2 = Y^r(5/2 + i)-^ 



V (^^ + >^ + f ^ + f )r(m + |)r(. + 1 + V2) _ 

^,„r(z + i)r(m + 1) ^3 I J 



where S'fcg = {k\ - !)/(! - /cs)^, S'i.m = + 2(i + m) + 5, and 7„ii„ is given by 
Eq.(34). 

The differences between the calculation with a simple electron energy distribu- 
tion /(7) = 7^"^ and the correct distribution 7(7) = are important 
for ei/e < 20 (Fig.l7). 
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